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Abstract 

We show how string breaking on the lattice, treated as a mixing effect, can be 
related to decay rates for heavy quark systems. We use this to make a preliminary 
calculation of the energy split at maximum mixing for static quarks in QCD from 
the decay rate for Y(4S) — ► BB . We extend the calculation to achieve rough 
estimates for the contributions of channels involving B, B* , B s and B* mesons 
to the width of the T(55) . 
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1 Introduction 



It has been recognised for a long time that, just as the confining string-like interquark 
potential provides a good basis for understanding valence quark dynamics, so the break- 
ing of the string provides a natural way of understanding the influence of sea quarks 
on decay processes jl], ||. Lattice measurements of the confining QCD potential, in 
quenched and unquenched calculations, are well understood [§, §, |[ §, [7|, f| f|. More- 
over the simplest type of quark model using this potential gives a reasonable account of 
the observed heavy meson spectrum. Recently there has been progress in the attempt 



to observe string breaking effects directly on the lattice [1C, 11, 12]. 
In this paper we investigate how sufficiently accurate lattice measurements of the in- 
terquark potential and string breaking effects might be combined to implement a pro- 
gramme of computation of decay rates for heavy quark systems. The analysis depends 
on a two-channel mixing picture of the string breaking process suggested by a crude 



model proposed previously [13], [14], |T3] . The basis of the mechanism is the fact that the 



light quark-anti-quark pair is produced (or absorbed) in a triplet spin state. Although 
our picture of string breaking is different in detail from that of Isgur and Paton [|I]] the 
results are still consistent with the phenomenological 3 Po- m odel . 

2 Crude Model 



The crude model described previously |T^, [14], |T§], was designed to provide a simple 



two-channel mixing model of string breaking. It suggests that when the string breaks 
each light quark and anti-quark materialises in the neighbourhood of the appropriate 
static anti-quark or quark. In this respect it differs from the string breaking model 
of Isgur and Paton jl] who viewed the materialisation of the light quark-anti-quark 
pair as a local event on the string. As a justification of our picture we point out 
that to a first approximation the local gauge field energy density does not change as 
the static quarks are separated. There is no real mechanism therefore for the local 
production of the light pair. The breaking event is only possible when the total energy 
distributed along the string is sufficient to support light pair creation. Our model then 
suggests that this creation occurs essentially instantaneously in such a way that the 
gauge field distribution supported by the static quarks is cancelled by the opposing 
gauge field structure of the light quarks. It is this idea of instantaneous transition that 
we exploit below in our development of the model. However the outcome of the these 
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considerations is not in contradiction with the 3 Po model that is supported by the 
analysis of Isgur and Paton, at least for the decay processes to which we apply our 
ideas. 

As was pointed out in the original description of the crude model O, O, the 



dynamical quark-anti-quark pair are produced in a correlated spin state. Here we 
examine this mechanism a little more closely. In the model the strong coupling graph 
that corresponds to the production of the pair is shown in Fig 1 . The quark or anti- 
quark binds in an S'-wave with the appropriate static anti-quark or quark. Choosing 
as the (imaginary) time direction and 1 as the direction of spatial separation of the 



static quarks, we see that the rules for quark propagation |Lj|, [14], [y|, [L9| gQj yield an 
amplitude with spin structure of the form 



1 + 7o 1 + 7i 1 - 7o _ 1 / (7i 
2 2 2 " 2 I 



This can be interpreted as meaning that the light quark-anti-quark pair are produced 
with a spin wave function 

£a/3 = -^{Vi^HjaP , (2) 

where a and (3 are the spin labels for the quark and anti-quark respectively and R 
is the unit vector pointing along the line of separation of the the static quarks. This 
shows that the light quark- ant i- quark pair are produced in a triplet state with zero 
spin component along the line of separation of the static quarks. 

3 Basic Mixing Scenario 

In order to explain the basic idea we first consider a model in which the particles of 
the theory have no spin but do experience confinement. The non-abelian Higgs models 



are of this type [TO, ITTJ . Moreover they exhibit string breaking on the lattice of a type 



consistent with a two-channel mixing scenario [|IB|, [TJ], 15]. The channels involved are 

1) a static quark-anti-quark {QQ} system with a separation R and with a connecting 
flux string 

2) a heavy-light meson-meson system {Qq, qQ} . 

For the moment, we ignore the spins of the quarks and treat them as scalar particles. 
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On the lattice it is possible to measure the string-string correlator QssiT) as the RxT 
Wilson loop. It has the form 

g ss {T) ~ e- v QQ^ T , (3) 

where Vqq(R) can be interpreted as the interquark potential experienced by the static 
quarks. Lattice measurents are consistent with the Cornell potential |T^, [H], [Ll| 

V QQ (R) = C-^ + aR , (4) 

where a is the string tension. The correlator for two-meson states over the time interval 
T can be obtained from measurements of observables appropriate to the two-meson 
state. Its behaviour is of the form 

W(T)~e- £ —^ T , (5) 

where E MM (R) = 2E M + Vmm{R) and E M is the energy of a single static meson and 
Vmm{R) is the potential interaction between the two static mesons. This potential en- 
ergy has not been measured in QCD but does show up in the Higgs model calculations. 



We might expect it to have a form like [15 



Vmm(R) = W M e- mR , (6) 

where m is the mass of a light meson that can be exchanged between the heavy mesons 

and Wm is an energy determining the overall strength of the interaction. 

This description is of course an oversimplification in the presence of dynamical matter 



fields. This is clear from the Higgs model calculations [TO, O] and in QCD from the 



thermal Polyakov loop calculations The reason is that string breaking can occur as 
a process that leads to a mixing of the string and two-meson states. Consistently with 
these static model calculations we may represent the mixing process by a transition 
potential Vi(R) . On the basis of the crude model we might expect it to have the form 

0, MM 

Vj(R) = We~ m " R , (7) 

where m q is the mass of the dynamical matter field and W is the energy parameter 
determining the strength of the transition potential. It follows that we are dealing with 
a two-channel interaction matrix of the form 
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The energies measured directly on the lattice are the eigenvalues of V(R) . They are 

V ± (R) = \ [V QQ {R) + E MM {R) =F y/(V QQ (R) - E M m(R)) 2 + 4(^(i?)) 2 } . (9) 

We keep the convention used previously that V+(R) is the lower of the two eigenenergies. 
If Vj(R) is sufficiently small then the eigenenergies are dominated by the diagonal 
elements of V(R) . We can define a critical value of R c for which these diagonal values 
are equal 

V QQ (R C ) = E MM (R C ) . (10) 
When R = R c the the split in the eigenenergies is proportional to Vj(R c ), we have 

V ± (R c ) = V qQ (R c )tVj(R c ) . (11) 

For R « R c we expect Vqq(R) << E M m{R) with the result V+(R) ~ Vqq(R) and 
V-(R) ~ Emm(R) ■ For R >> R c we expect V q q(R) » E MM (R) with the result 
V-(R) — Vqq(R) and V + (R) ~ Emm{R) ■ The interchange that occurs near R = R c 



is the mixing phenomenon. Just as in refs |13], |14[ y| we can compute a mixing angle 
6* given by 



-(V Q(i2) - Emm {R) ) + v / (^q(^) " Emm(R)) 2 + ^{R)? 

taXLV — ; : . 12 

2Vj(R) V ; 

We have 

^=0(^^)0-, (13) 

where 

„ / cos^ — sin^ \ 
V smu cosy / 

As i? passes from below to above R c , the mixing angle increases from to tt/2 hitting 
7r/4 when R = R c . The goal then is to use lattice measurements to determine the 
eigenenergies and the independence of the mixing angle. By those means we can 
determine the full two-channel interaction matrix, V(R) . We discuss this in more 
detail below. 

Once we know the elements of V(R) we can make the two-channel equivalent of 
the Born-Oppenheimer approximation and replace the static quarks by heavy non- 
relativistic quarks and mesons. The effective Schrodinger equation is 

z^V(R) = ( ~^J V2 _ J_ v2 ) V(R) + V(R)iP(R) , (15) 
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where 

^ = ( ) ' < 16 » 

and hq and \xm are the reduced masses for the heavy quarks and heavy mesons in their 
respective channels. These also can be obtained from appropriate lattice measurements 
of heavy single meson states. 

We can exploit eq fll5f ) by first ignoring the off-diagonal elements of V(R) . We then 
calculate the heavy QQ bound states with energies Eqq . If this energy lies above the 
two-meson channel threshold then we can compute the decay rate as 

r = 2np(k)\T(k)\ 2 , (17) 

where Eqq is the energy of the bound state QQ system in the potential Vqq(R) , k is 
the relative momentum of the final state mesons, p(k) is the density of states factor. 
The energy of the two-meson state is E(k) which satisfies 

E(k) = -^-k 2 + 2E M = E Q Q- . (18) 

The transition amplitude is 

T(k) = (k\Vj\^ Q ) , (19) 

|k) being the wavefunction, with incoming scattered wave, for the two meson system in 
the presence of the potential Vmm(R) • As a first approximation we could neglect the 
the final state interactions represented by the scattered part of the the wave function 
and replace |k) with the plane wave e lk R . We have then 

T(k) = J ^Re-^V^MR) = ^J o °°dRR 2 Vj(R)4j q (R) . (20) 

In the above it is assumed that the static approximation for the mesons is accurate in 
the sense that 

Eqq-2E m = M Q q-2M m . (21) 

where Mqq and Mm are the actual masses of the heavy quark state and the heavy-light 
meson state respectively. To the extent that the results are only approximate the left 
side of eq([21|) should be replaced by the right side in an actual estimate of the decay 
rate. 
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4 Lattice Measurement of the Interaction Matrix 



On the lattice we can measure V±(R) [It], [II], |T2| . If we are confident of our parametriza- 
tions of the elements of V(R), then it is likely that we can determine them by perform- 
ing a fit to the eigenenergies provided they are known sufficiently accurately and on a 
sufficiently fine spatial lattice. 

In principle, the the eigenenergies are not sufficient to fix the elements of V(R) . A 
superior approach, therefore, would be to measure also the mixing angle 9 . The results 
obtained in the S£/(2)-Higgs models for the independence of the ground state overlap 



of the lower eigenmode suggest that this may well be possible ||10| , |Tl| . The two-channel 
model suggests that the contribution of the lower energy eigenmode to the Wilson loop 
is 

Gss(R) - Tcos6e~ v+(R)T Tcos6 , (22) 

where T is independent of R . If we make this assumption then then we can deduce 
the independence of 9 from the measurement of Qss(R) ■ 

It would be interesting if this scenario could be generalized to a set of observables 
Ai(R, T) i = 1, . . . , N for which the correlators have the form 

(MR,T)MR,0)) = £ T ta (R)e- Va(R)T T ja (R) , (23) 

where the independence of the coupling coefficients comes from the movement of the 
mixing angle. That is 

T lQ (R)= E ^iaO^R) • (24) 

a=S,M 

Introducing the iV-component vectors C a = {T ia } and C a = {T ia } we easily see from 
the the form of the orthogonal matrix O given in eq flli]) that 

C + ( J R)AC_( J R) = C 5 AC M , (25) 

and 

C l+ {Rf + CUR? = C 2 lS + C 2 m i = l,...,JV . (26) 

If the left sides of these equations do show a lack of dependence on R then we would 
have a good test of the the above coupling hypothesis. It is not obvious that all 
sets of operators have this property but it may be possible to construct a sufficient 
set. However we might reasonably expect in general that the i?-dependence of these 
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quantities is relatively weak compared to that induced by the rapid variation of the 
mixing angle through the mixing region. 

On the assumption that an acceptable set of operators can be found, the dependence of 
9 on R can be elucidated by examining the quantities y + = n.C + (R) and y_ = n.C-(R) 
where n is any appropriate iV-vector. They have the form 

V+ \ = ( cos0 -sin 6> \ / y s 
V- J V sin6> cos6> J \ y M 

where ys = n.Cs and yM = h.Cm • Clearly the point (y + ,y_) lies on a circle and 9 is 
the angle, referred to an appropriate origin, that fixes the position of the point on the 
circle. 




5 Light Quarks with Spin 

Static quark calculations are particularly relevant to the spectator quark approximation 
in which, except for determining multiplicities, the heavy quark spins play essentially 
no role. Of course in the physical case of bottomonium, the small mass differences 
associated with the hyper-fine structure of the 5-mesons are crucial in determining 
allowed decays. These effects are however essentially kinematic and will be taken 
into account at the appropriate point. In examining the basic phenomenon of string 
breaking we will temporarily omit these spins from consideration. The light quark 
spins however play a crucial part in both string breaking and decay dynamics. 
As discussed in section 2 our crude model provides a generalization of the scalar case 
to that of light quarks with spin. It remains a two-channel mixing problem because 
the light quark-anti-quark pair emerges in a triplet spin state with definite orientation. 
We can achieve this result by introducing a transition potential of the form 

™ - ( as ) ■ «> 

where now the spin label structure of the transition element is given by 

V MSaP (R) = ^VKi?)(<7iRi)a/3 • (29) 

This indeed implies that Vms(R) couples the static quark- ant i- quark state to a two- 
meson state in which the light quark spins are in a state represented by the wave 
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function £ al 3 given in eq(^) . Projecting down onto the appropriate subspace the two 
channel structure emerges and the effective potential matrix is 

This is of precisely the form of the interaction matrix in the scalar case and the eigen- 
values V±(R) which which are given by eq(|9]) are the energies that will be computed on 
the lattice from the measurement of appropriate operators. In the same spirit as for the 
scalar calculation we can expect, through careful lattice measurements, to determine 
the entries in V (R) . Having done that we revert to considering the Schrodinger 
equation 

i^(R) = ( ) V(R) + V(K)^(R) , (31) 

\ 2/j.m / 

Following the pattern set out for the scalar case we calculate the decay rate by solving 
eq fl3"ID with Vi(R) = for an S'-wave bound by Vqq(R) and a scattering wave controlled 
by Vmm(R) ■ The transition rate is then 

r = 2vr£p(A;)|7Xk)| 2 , (32) 

i 

where 

T l (k) = (ki\V MS (R)\^ Q ) , (33) 

and i is the polarization state of the triplet spin wavefunction. In principle the state 
|k i) should be a solution of the the scattering problem with the potential Vmm(R)- If 
we make the approximation of neglecting final state interactions we can set 

|k i) = -^cr^ . (34) 

We have then 

Ti(k) = -j= J ^Re-^TricXiVMsimMR) ■ (35) 



s/2 

That is 

T<(k) = 47^*, f dRR (°-^ - Vj{R)^ Q {R) . (36) 



Note that Tj(k) ~ 0(k) for small k . This is the angular momentum barrier associated 
with the P-wave final state. 
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6 Heavy Meson Decay 



Currently there is no detailed information on string breaking in QCD at zero tempera- 



ture although there is a strong indication of the mixing effect at finite temperature fL2 |. 
The latest measurement of the QQ potential with dynamical quarks does not show an 
unequivocal mixing effect directly ||. This is not inconsistent with the mixing scenario 
if the mixing range in R and the split in eigenenergies at maximal mixing is sufficiently 



small [0, [14|, . It is therefore interesting to reverse the above argument and use the 



information on the decay of the T(4S) to BB to estimate the energy split at maximum 
mixing that we should expect to see on the lattice. The theory also enables us to make 
predictions for other bottomonium decays. 

The hyperfine structure of the heavy-heavy and heavy-light quark systems has a strong 
effect on the relative positions of masses and thresholds and hence on decay rates. In 
order to make use of experimental results therefore, it is necessary to reinstate the 
spins of the heavy quarks. To compute the decay rates of real processes, we must take 
into account the spin structure of the open final state channels in order to assign the 
correct momenta to them. 

The heavy quarks in the T are in a triplet state and continue, as spectator quarks, in 
that state after decay. We can represent this state by a spin wave function crja'p' /v2, 
where a' is the heavy quark label and f3' is the heavy anti-quark label. In our model the 
light quarks are produced in a triplet state that can be represented by a wave function 
Piap/Vz ■ A sufficiently complete basis for the achievable final states is 

^ij = 2 a i<*l3 (T joi>l3' ■ (37) 

The transition matrix element we require is then 

(ipif f \s.\V M s\ipQ,j} = SjfTif^k) , (38) 

where k is the relative momentum for the final state mesons and Ty (k) is the amplitude 
as previously defined. The spin state ipij can be decomposed into a physical basis such 
as BB, BB* + B*B, B*B*(S = 2) and B*B*(S = 0) . The resulting decay rate is 
given by eq (|32|) modified by a probability factor p appropriate to the channel under 
consideration. 

T = p27rp(k)Y / \T l (k)\ 2 , (39) 

i 

These probabilities, calculated from the appropriate overlap coefficients in the spin 
wavefunction recoupling scheme, for the physical channels are 
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BB 


BB* + B*B 


B*B*(S = 2) 


B*B*(S = 0) 


1/12 


1/3 


5/9 


1/36 



Introducing the radial wavefunction x(R) = v4^R4>q{R) and using eq(0) we can 
compute Tj(k) from eq(|36]) as 



Ti(k) = V47rifciA(fc,m 9 )W , (40) 

where 

am \ f°° 7T ^f cos kR sin kR\ _ m R . 

m s ) = / o dii ^— - — j a -^(fl) . (41) 

Finally the decay rate to a particular channel is 

r =pN f 4fik 3 \A(k,m q )\ 2 W 2 J (42) 

where we have used p(k) = \ikj (27r 2 ), // is the reduced mass of the final state particles 
and Nf is the number of degenerate BB channels which is identical to the number of 
degenerate light quarks. We use mass values from the Review of Particle Physics [pT] , 



Mass 


M B 


M B * 


M Bs 


Mb. 


Mr (as) 


Mr(5S) 


GeV 


5.279 


5.325 


5.369 


5.416 


10.580 


10.860 



We choose as a representation of Vqq(R) its parametrization in eq(|). A choice of 
parameter values that accord reasonably well with lattice measurements and the T(S) 
excitation spectrum is a = 0.52 and o = (0.429) 2 GeV 2 , together with a heavy 
quark reduced mass //q = 2.1 GeV. The resulting radial wavefunctions x(R) f° r T(4S r ) 
and T(6S) are shown in Fig 2. Both wavefunctions remain substantial over a range of 
10 GeV -1 ~ 2 fm with oscillations near the origin of a half wavelength ~ 1 — 2 GeV -1 
before reaching broad maxima at ~ 5 & 7 GeV -1 . 

As indicated above we ignore final state interactions, effectively setting Vmm(R) = . 
Because our calculation is a preliminary one we will assume in what follows that m q = 
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. This may be justified as a rough approximation by the low masses of dynamical 
quarks. A crucial factor in computing the decay rate is the amplitude A(k,m q ). We 
exhibit this amplitude, for the case m q = 0, as a function of the final state momentum 
k in Fig 3 for the T(4S) and T(5S). Note that in both cases it is small in magnitude 
for k > 0.5 GeV. We can expect therefore as a general rule that only decays with a 
sufficiently small final state momentum will contribute substantially to the width of 
the initial state. 

In our model the only open channel for T(4S')-decay is BB . Experimentally this 
process dominates. The relative momentum of the the BB mesons in the final state is 



k = V /M B (M 4S - 2M B ) = .3408GeV . This yields a value for the amplitude A(k, 0) = 
-10.4 . From eqflU) with N f = 2 and p = 1/12 we find 

T(T -> BB) = 7.02 W 2 . (43) 

The experimental result is T(T — > BB) = 0.010(5) GeV. This yields an estimate 
W = 0.038(9) Gev. In turn we can estimate the eigenenergy split at maximum mixing 
as AE = 2^J NfW = 0.11(3) GeV . We can crudely estimate the range in R over which 
the mixing takes place as 

AE AE 

with the result AR ~ 0.6 GeV -1 ~ 0.12fm . To resolve such a split requires a rather fine 
lattice. The small value of this split, coupled with the decoupling of the two-meson 
state outside the mixing region may be regarded as an explanation of why current 
measurements of the Wilson loop do not yet reveal string breaking in lattice QCD. 
We can extend the above analysis to higher excitations such as the T(5S*) which has 
a number of open decay channels. The results, assuming W has the same value as 
before, are 
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Channel 


N, 


k(GeV) 


A(k, 0) 


r(MeV) 


BB 


2 


1.26 


-0.68 


2.0 


BB* + B*B 


2 


1.16 


-0.734 


9.0 


B*B* 


2 


1.06 


-0.430 


4.0 


B S B S 


1 


0.81 


1.651 


2.0 


B S B* S + B* S B S 


1 


0.636 


1.1545 


3.0 


B* S B* S 


1 


0.389 


-7.33 


29.0 



The total width is T ~ 50 MeV which is to be compared with the the quoted ex- 
perimental result T = 110(13) MeV. Given the crudity of the theoretical approach and 
the experimental difficulties this is not an unreasonable comparison. 

7 Conclusions 

We have examined the relevance of string breaking to the decay of heavy quark systems. 
The anaysis shows how string breaking on the lattice can be related to the decay rates 
for such processes as T(AS) — > BB . The central idea is that the string breaking can be 
viewed as a two channel process. This is controlled by a 2 x 2 potential matrix of which 
the diagonal elements represent the interquark potential and the energy and interaction 
potential of two static mesons. The off-diagonal elements correspond to a transition 
potential between the two channels. Just as the interquark potential can be used in a 
dynamical calculation of heavy quark states so the transition potential can be used to 
compute the transition rate from the bound quark system to two freely moving heavy 
mesons. This is a standard relatively non-controversial quantum mechanics calculation. 
The main assumption behind it is that the Born-Oppenheimer approximation that 
underpins the heavy quark calculation holds good when the transition matrix elements 
are included. This is equivalent to assuming that there is no strong retardation effect 
in the transmission of the dynamical quark. If this quark mass is sufficiently low this 
assumption is not unreasonable. 

At the moment there are no measurements in QCD that permit an accurate realization 
of the scheme. However, with some simplifying assumptions on the form of the elements 
of the potential matrix, it is possible to make a preliminary calculation of the decay rate 
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of the T(4S) in terms of the energy parameter W, yielding a value ~ 38 MeV. This gives 
rise to an eigenenergy split of ~ 110 MeV. Although this suggests that the split may 
be hard to detect on the lattice the study of mixing remains an important goal because 
of its physical significance. Our analysis also strongly supports the desirability of 
measuring an appropriate suite of operators for the purpose of measuring the movement 
of the mixing angle fl(| 11 1. In fact, as indicated in [15|, it would be very useful in 



estimating Vj(i2) to measure the transition amplitude for string to two mesons in the 
quenched appoximation ||22|| . 

Taking the above value of W as a guide we are able to use the formalism to compute 
the partial widths for T(5S)-decay. We find that the dominant decay is through the 
B*B* channel. Our computed total width is roughly half the measured width. We feel 
that this result is not unreasonable given the preliminary character of the calculation 
and the limited information available from experiment on decay channels. 
We conclude that with some improvement of lattice measurements it will be possible 
to confront string breaking with experimental results when these become available. 
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Figure 1: Transition from string to two-meson state. Heavy lines represent static 
quarks, light lines represent dynamical quarks. 
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Y: 4S (solid) 5S (dashed) 



DC 




R (GeV) 



Figure 2: Wavefunctions for T(AS) and T(5S). 
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Figure 3: Transition matrix element, A(k,0), for T(AS) and T(5S). 
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